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Abstract. The enumeration of perfect matchings of graphs is equivalent to 
the dimer problem which has applications in statistical physics. A graph G 
is said to be n-rotation symmetric if the cyclic group of order n is a sub- 
group of the automorphism group of G. Jockusch (Perfect matchings and per- 
fect squares, J. Combin. Theory Ser. A, 67(1994), 100-115) and Kuperberg 
(An exploration of the permanent-determinant method. Electron. J. Combin., 
5(1998), #46) proved independently that if G is a plane bipartite graph of 
order A'^ with 2n-rotation symmetry, then the number of perfect matchings of 
G can be expressed as the product of n determinants of order N/2n. In this 
paper we give this result a new presentation. We use this result to compute 
the entropy of a bulk plane bipartite lattice with 2n-notation symmetry. We 
obtain explicit expressions for the numbers of perfect matchings and entropies 
for two types of cylinders. Using the results on the entropy of the torus ob- 
tained by Kenyon, Okounkov, and Sheffield (Dimers and amoebae, Ann. Math. 
163(2006), 1019-1056) and by Salinas and Nagle (Theory of the phase tran- 
sition in the layered hydrogen-bonded SnCfi ■ 2H2O crystal, Phys. Rev. B, 
9(1974), 4920-4931), we show that each of the cylinders considered and its 
corresponding torus have the same entropy. Finally, we pose some problems. 



1. Introduction 

An automorphism of a graph G is a graph isomorphism with itself, i.e., a bijection 
from the vertex set of G to itself such that the adjacency relation is preserved. The 
set of automorphisms defines a permutation group on the vertices of G, which is 
called the automorphism group of the graph G. The automorphism group of a 
graph G characterizes its symmetries, and is therefore very useful for simplifying 
the computation of some of its invariants. One of the most successful results in 
this direction is to use the character of the automorphism group of a graph G to 
compute its characteristic polynomial (see for example Cvetkovic et al. [5|). Recall 
that a cyclic group is a group that can be generated by a single element x (the 
group generator). A cyclic group of order n is denoted by C„, and its generator x 
satisfies = 1, where 1 is the identity element. A graph G is said to be n-rotation 
symmetric if the cyclic group C„ is a subgroup of the automorphism group of G. 
We call Cn the rotational subgroup of an n-rotation symmetric graph G. We say 
that a plane graph G is reflective symmetric if it is invariant under the reflection 
across some straight line. 
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A perfect matching of a graph G is a cohection of vertex-disjoint edges that 
are collectively incident to all vertices of G. Let M{G) denote the number of per- 
fect matchings of a simple graph G. Problems involving enumeration of perfect 
matchings have been examined extensively not only by mathematicians (see for 
example [2l[3l|4l[8l[T4l[23l[26]) but also by physicists and chemists (see for example 
[SlIIIlinilialEllES]). In 1961, Kasteleyn TI] foimd a formula for the number of 
perfect matchings of an m x n quadratic lattice graph. Temperley and Fisher [27] 
used a different method and arrived at the same result at almost exactly the same 
time. Both lines of calculation showed that the logarithm of the number of perfect 
matchings, divided by converges to 2c/7r k, 0.5831 as m,n ^ oo, where c is 
Catalan's constant. This limit is called the entropy of the quadratic lattice graph 
and the corresponding problem was called the dimer problem by the statistical 
physicists. In 1992, Elkies et al. 6J studied the enumeration of perfect matchings 
of regions called Aztec diamonds, and showed that the entropy equals « 0.35. 
Cohn, Kenyon, and Propp demonstrated that the behavior of random perfect 
matchings of large regions R was determined by a variational (or entropy maxi- 
mization) principle, as was conjectured in Section 8 of 6J, and they gave an exact 
formula for the entropy of simply-connected regions of arbitrary shape. Particularly, 
they showed that computation of the entropy is intimately linked with an under- 
standing of long-range variations in the local statistics of random domino tilings. 
Kenyon, Okounkov, and Sheffield fl4J considered the problem enumerating perfect 
matchings of the doubly-period bipartite graph on a torus, which generalized the 
results in [4 . They proved that the number of perfect matchings of the doubly- 
period plane bipartite graph G can be expressed in terms of four determinants and 
they expressed the entropy of G as a double integral. 

The dimer problem on hexagonal lattice graphs has also been examined in the 
past [3 [31] . It can be seen as equivalent to the combinatorial problem of "plane 
partitions" (see MacMahon [IS])- The a, 6, c semiregular hexagon is the hexagon 
whose side lengths are in cyclic order, a, 5, c, a, 6, c. Lozenge tilings of this region 
are in correspondence with plane partitions with at most a rows, at most h columns, 
and no part exceeding c. MacMahon 19J showed that the number of such plane 
partitions (i.e., the number of perfect matchings of the corresponding hexagonal 
lattice graph G(a, 6, c)) equals 

lilili,+,+,_2 

Elser [7] obtained the entropy of G{aTa,a) which is approximately 0.2616. The 
dimer problem on some other types of plane lattices are also considered (see Wu 
[30] and the references cited therein). 

Jockusch [TU] first used a combinatorial method to obtained a product formula for 
the number of perfect matchings of a type of rotational symmetric bipartite graphs. 
Kuperberg [T5] used the representation theory of groups to obtain independently a 
product formula which can be equivalent to the Jockusch's one. Ciucu [2] , and Yan 
and Zhang [32l [33| [34] considered the enumeration of perfect matchings of general 
graphs (lattices) with a certain type of reflective symmetry. Ciucu [2] obtained a 
basic factorization theorem for the number of perfect matchings of plane bipartite 
graphs with reflective symmetry. Yan and Zhang [32l [33] extended Ciucu's result 
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to general plane graphs (not necessary to be bipartite) with reflective symmetry by 
using PfafRans. 

The current paper deals with the enumerating problem for perfect matchings of 
cylinders which can be regarded as the graphs with rotational symmetry. In the 
next section, we introduce the Pfaffian method for enumerating perfect matchings 
of plane graphs (see for example [3 [HI [26l El])- In Section 3, we give a new 
presentation of the product formula for the number of perfect matchings of plane 
bipartite graphs with 2n-rotation symmetry, which was found independently by 
Jockusch [10 and by Kuperberg [T5^. We use this result to compute the entropy of 
a bulk plane bipartite lattice with 2n-notation symmetry. In Section 4, we obtain 
the explicit expressions for the numbers of perfect matchings and entropies for two 
types of cylinders (bipartite graphs with the cylinder-boundary condition). Using 
the theorem obtained by Kenyon, Okounkov, and Sheffield [M], we compute the 
entropy of a type of toruses (bipartite graphs with the torus-boundary condition) . 
On the other hand, the entropy of another type of toruses was computed by Salinas 
and Nagle [25]. Our results show that each of the cylinders considered and its 
corresponding torus have the same entropy (a similar result for the plane quadratic 
lattices had been obtained by Kasteleyn [H]). Finally, Section 5 poses two open 
problems. 

2. Pfaffians 

Let B = (bij)2nx2n bc a skew symmetric matrix of order 2n. For each partition 
P = {{h,ji},{i'2,j2},---,{'injn}} of the set {l,2,...,2n} into pairs, form the 
expression 

bp = sgn{iijii2j2 ■ ■ ■ injn)bitjtbi^j2 ■ ■ ■ bi„j„, 

where sgn(iijii2j2 ■ ■ ■ injn) denotes the sign of the permutation iijii2j2 ■ • ■ injn- 
Note that bp depends neither on the order in which the classes of the partition are 
listed nor on the order of the two elements of a class. So bp indeed depends only on 
the choice of the partition P. The Pfaffian of the skew matrix B (see [T7|), denoted 
by Pf{B), is defined as 

Pf{B) = J2bp, 
p 

where the summation ranges over all partitions of {1, 2, . . . , 2ri} which have the 
form of P. 

Theorem 2.1 (The Cayley's Theorem, fTT]). Let B = (6ij)2nx2ri be a skew sym- 
metric matrix of order of2n. Then 

det(B)^Pf{Bf. 

The Pfaffian method for enumerating perfect matchings of plane graphs was in- 
dependently discovered by Fisher , Kasteleyn [11] , and Temperley [27] . Given a 
plane graph G, the method produces a matrix A such that the number of perfect 
matchings of G can be expressed by the determinant of matrix A. By using this 
method, Fisher [8j, Kasteleyn [11], and Temperley [27] solved independently a fa- 
mous problem on enumerating perfect matchings of an m x n quadratic lattice graph 
in statistical physics-Dimer problem. Given a simple graph G = {V{G), E(G)) with 
vertex set V{G) — {ui,W2, . . • ,wjv}, let be an arbitrary orientation. The skew 



4 



WEIGEN YAN, YEONG-NAN YEH, AND FUJI ZHANG 



adjacency matrix of G"^, denoted by A{G'^), is defined as follows: 
where 

{1 if {vi,Vj) is an arc of 
— 1 if {vj, Vi) is an arc of C, 
otherwise. 

Obvious, A{G'^) is a skew symmetric matrix. 

If D is an orientation of a graph G and C is a cycle of even length, we say that 
G is oddly oriented in I? if C contains odd number of edges that are directed in D 
in the direction of each orientation of G. We say that D is a Pfaffian orientation 
of G if every nice cycle of even length of G is oddly oriented in D (a cycle C of G 
is called be nice if the induced subgraph G — G of G by V{G)\V{G) has perfect 
matchings). It is well known that if a graph G contains no subdivision of K^^^ 
then G has a Pfafiian orientation (see Little [E]). McCuaig [20], and McCuaig, 
Robertson et al [21], and Robertson, Seymour et al f2T found a polynomial-time 
algorithm to show whether a bipartite graph has a Pfaffian orientation. Stembridge 
|26j proved that the number (or generating function) of nonintersecting r-tuples of 
paths from a set of r vertices to a specified region in an acyclic digraph D can, 
under favorable circumstances, be expressed as a Pfaffian. See a survey of Pfafffian 
orientations of graphs in Thomas |28j . 

Lemma 2.2 (Lovasz et al. [I? ). //C^ is a Pfaffian orientation of a graph G, then 

M{G) = \Pf{A{G'))\, 
where A(G'^) is the skew adjacency matrix of G'^ . 

Lemma 2.3 (Lovasz et al. [17j). If G is a plane graph and G*^, an orientation of G 
such that every boundary face-except possibly the infinite face-has an odd number of 
edges oriented clockwise, then in every cycle the number of edges oriented clockwise 
is of opposite parity to number of vertices of G'^ inside the cycle. Consequently, G'^ 
is a Pfaffian orientation. 

Lemma 2.4 (J8\ llll [27] ). Every plane graph has a Pfaffian orientation satisfying 
the condition in Lemma 2.3. 

3. A PRODUCT THEOREM 

Let G = {V{G),E{G)) be a simple connected plane bipartite graph of order 
N with 2n-rotation symmetry, which has symmetric axes Iq, £i, . . ., and ^2n-i 
passing through the rotation center O. Figure 1(a) shows an example of a 6- 
rotation symmetric graph. If G has some vertices lying on its symmetric axes ^i's, 
then we can rearrange vertices of G such that there are no vertex lying on its new 
symmetric axes. Hence we can assume that there exists no vertex of G lying on Iq, 
£i, . . ., ^2n-i- If we delete all the edges intersected by £q, £i, . . ., and £2n-i, then 
2n isomorphic components Gi = {V{Gi), E{Gi)) for < i < 2n — 1 are obtained, 
where V{Gi) — {vi^\v2\ . ■ . ,v'"j^^^^} is the vertex set of Gi. Note that G is a 
plane graph. Hence there exists no edges between Gi and Gj for j ^ i — l,i + 1 
(mod 2n). That is, we can arrange all G^'s on a circle clockwise (see Figure 1(a)). 
Without loss of generality, we assume that the set of edges of G between Gi and 



MATCHINGS AND ENTROPIES OF CYLINDERS 



5 
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(a) (b) 

Figure 1. (a) A connected plane bipartite graph G with 6- 
rotation symmetry, where the local width w{G) of G equals 2. 
(b) One component Gq of G. 



Gi+i is exactly {vrjvi[ '\k = l,2,...,p} (see Figure 1(a)), where G2n = Go, 
Wsfc"'' = vi^J , and Wrf"'' = vi'f^ , and the generator x of the rotation subgroup C2n 
of the automorphism group of G maps Vr^ (resp. ui*'') to Vr^^^ (resp. Vs^^^^) 
for I < k < p, respectively. The local width of a 2n-rotation symmetric graph G 
with 2n isomorphic components G^'s, denoted w{G), is defined to be the number of 
edges in Gi lying on the boundary face of G which contains the rotation center O. 
For the graph G illustrated in Figure 1(a), w{G) — 2. For the underlying graph G 
(which is a 6-rotation symmetric graph) of the digraph illustrated in Figure 4(a), 
w{G) = 1. 

Note that Go is a plane graph. By Lemma 2.4, there exists a Pfafhan orientation 
Gq of Go satisfying the condition in Lemma 2.3. For the corresponding graph Go 
illustrated in Figure 1(b), the orientation Go shown in Figure 2(a) is a Pfaffian 
orientation satisfying the condition in Lemma 2.3. Now we define a Pfaffian orien- 
tation G*^ from Gq as follows. If we reverse the orientation of each arc of Gg, then 
we can obtain an orientation of Go satisfying conditions in Lemma 2.3, denoted by 
Gq'^. Hence, Gq'^ is the converse of Gq. Note that Go is bipartite. Therefore, Gq"^ 
is a Pfaffian orientation of Go satisfying the condition in Lemma 2.3. Since all G^'s 
are isomorphic, we define the orientations of all G2i's to be the same as Gg and the 
orientations of all G2i+i's to be same as Gq"^ for < i < n — 1, respectively. In 
order to obtain an orientation G^ of G, we need to orient all edges with the form 
of vr^v^^'^^^ in G. Note that the subgraph G04 of G induced by V{Go) U V{Gi) 
is a plane graph. Define the orientations of Gq and Gi in Go,i to be Gg and Gq"^, 
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(a) (b) 

Figure 2. (a) A PfafRan orientation Gq of Go- (b) The corre- 
sponding Pfaffian orientation Gg ^ of Go,i. 



respectively, and orient edges between Gq and Gi one by one such that every face 
between Go and Gi has an odd number of edges oriented clockwise. Hence we have 
obtained an orientation Gq i of Go,i satisfying the condition in Lemma 2.3 such 
that the the induced orientations of Go and Gi in Gq i are exactly Gg and Gq"^. 
For the graph G illustrated in Figure 1(a) and the orientation Gq shown in Figure 
2(a), the corresponding orientation Gg ^ is pictured in Figure 2(b). 

Hence we have oriented the edges vi^Jvil'^ for 1 < fc < p in G. Define the 
orientations of edges vi^^vi^^^'' for 1 < i < 2n — 2 to be from vi^^ to vi'^^'^^ if the 
orientation of vi^Jvi]} is from vi^J to vi]} , and the orientations of edges fsl"'"^'' 
to be from vji'^^^ to vi^^ otherwise. So far, we have oriented all edges of G except 
the edges vl!l^~^'' vi^J for 1 < fc < p. In order to orient the edges vi!l"'~^\i^J such 
that G'^ is a Pfaffian orientation, we must distinguish the following four cases: 

Case (i). both n and w{G) (the local width of G) are even; 

Case (a), n is odd and u'(G) is even; 

Case (Hi), n is even and w{G) is odd; 

Case (iv). both n and w(G) are odd. 

For Cases (i), (ii), and (Hi), define the orientations of edges Wr^" ^^^s!} to be 
from vi^J^ to vi'f^ if the orientation of vi^Jvi]} is from vi]} to vi^J , and the orien- 
tations of edges vi^J^^^^'fJ to be from vi^J to w^^""^^ otherwise (see Figure 3(a)). 

For Case (iv), define the orientations of edges Wr^" ^"^vi^J to be from Wr^" to 
vi^J if the orientation of vi^^^^vi]^^ is from vi^J to vi]} , and the orientations of edges 
vi'l^^^^'vi^J to be from vi^J to fi-^""^^ otherwise (see Figure 4(a)). 

Hence we have obtained an orientation G^ of G. For the graph G in Figure 
1(a) and the orientation Gq of Go in Figure 2(a), by the above definition, the 
corresponding orientation G^ of G is shown in Figure 3(a). For the underlying 
graph G of the digraph illustrated in Figure 4(a), the corresponding orientation G^ 
is shown in Figure 4(a). Note that G is a bipartite graph. Every boundary face of 
the 2n Gj's in G"^ is oddly oriented clockwise. Moreover, it is not difficult to see 
that every boundary face in G^ intersected by the rotation axes £i,£2, ■ ■ ■ A2n-i- 
except the face containing the rotation center O-is oddly oriented clockwise. For 
Cases (i), (ii), and (Hi), the number of edges oddly oriented clockwise in the face 
containing the rotation center O equals w{G) x n + 2n — 1 = 1 (mod 2) (For the 
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Figure 4. (a) A regular orientation of G with the local width 
'w{G) = 1. (b) The corresponding weighted digraph Dj. 
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orientation illustrated in Figure 3(a), the number of edges oddly oriented clockwise 
in the face containing the rotation center O equals 2x3+5 = 11). For Case (iv), the 
number of edges oddly oriented clockwise in the face containing the rotation center 
O equals w{G) xn + 2n = 1 (mod 2) (For the orientation shown in Figure 4(a), the 
number of edges oddly oriented clockwise in the face containing the rotation center 
O equals 1x3 + 6 = 9). Thus the orientation of G satisfies the condition in 
Lemma 2.3. Consequently, G^ is a Pfaffian orientation of G. We call G^ a regular 
orientation of G. 

For j = 0,1, . . . ,n — 1, define aj (G) as 



a,(G) = 



cos ^ + i sin ^ i/ n = 1 {mod 2), w{G) = {mod 2); 

cos '"^''2n ^^ + * ^"^^tn^ Otherwise. 



(1) 

For convenience, set Uj ~ OLj{G) for < < n — 1. 

Now we construct n weighted digraphs Dq, Di, . . . , Dn-i of order N/2n from 
G^, which may contain both loops and multiple arcs, as follows. Note that Gg is 
a digraph containing neither loops nor multiple arcs and the vertex set of Gg is 
V{Gq) = U2°\ ■ • ■ , 'y^/2n}- ^0 the weighted digraph obtained from 

Gq by replacing each arc {vf'' , v^p ) with two arcs {vf^ , v^p ) and (uj"^ , vf'' ) with 
weights 1 and —1, respectively. For each j = 0, 1, 2, . . . , n — 1, define Dj to be the 

weighted digraph obtained from Gg by the following procedures: 

(1) For each i = 1,2, ... ,p, if wi"-* = v'f} , then add a loop at vertex ui^-* with 
weight aj+a~^ if (wi-^ '' , t'sj ) is an arc in G"^ and with weight —aj — a^^ if {v^s'} ,v'"r^) 
is an arc in G^. 

(2) For each i = l,2,...,p, if wi"' ^ vif' , then add two arcs (wr^'' . Ws^'' ) and 
{vi^\ Wr?'') with weights aj and aj^ if (wr?\ WsP) is an arc in G*^ and with weights 

—aj and —ctj^ if {vi]\vi.f^) is an arc in G*^. 

For the graph G illustrated in Figure 1 (a) , the corresponding weighted digraph 
Dj is shown in Figure 3(b), where each arc {vi^\vl.^^) of Dj, which is also be an arc 
in Gq , must be regarded as two arcs {vi^^ , v^^^ ) and {v^^^ , vi°^ ) with weights 1 and 
— 1, respectively. For the orientation G^ shown in Figure 4(a), the corresponding 
digraph Dj is illustrated in Figure 4(b). 

By a suitable labelling of vertices of the regular orientation G^ of G, for Cases 
(i)-(iii), the skew adjacency matrix Ai{G^) of G^ has the following form: 



^i(G^) 



A 


R 


•• 





rt 




-A 


R ■■ 











-R^ 


A •• 














•• 


A 


R 


-R 





•• 


■ -R^ 


-A 



2nx2n 
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and for Case (iv), the skew adjacency matrix ^2(6"^) of has the following form: 
/ 



V 



A 


R 


•• 





-R^ 




-A 


R ■■ 











~R^ 


A ■■ 














•• 


A 


R 


R 





•• 


■ -R^ 


-A 



where A is the skew adjacency matrix of Gq, and R is the adjacent relation in Gg i 
between Go and Gi. 

By the definitions of a^'s and Dj's, it is not difficult to see that the adjacency 
matrix Aj of weighted digraph Dj defined above equals exactly A + uj^ R + lo^^ R^ 



if n l(mod 2) and w(G) is even, and A + LOjR 



^R^ otherwise, where 



cos ^ + i sin ^ and 



(2j + l)ir , . . (2j + l)7r 



Now we can give a new presentation of the product theorem obtained indepen- 
dently by Jockusch jTfJ and by Kuperberg [Kj as follows. 

Theorem 3.1 (Product Theorem). Let G be a simple connected plane bipartite 
graph of order N with 2n-rotation symmetry and Gi 's for < i < 2n — 1 be the 2n 
graphs defined above. Suppose there are w(G) edges in Gi lying on the boundary 
face of G which contains the rotation center O and there exists no vertex of G lying 
on the rotation axes lo, h, ■ ■ hn-i passing O. Let G'^ be a regular orientation 
and Dj 's be the n weighted digraphs defined above. Then the number of perfect 
matchings of G can be expressed by 



M{G) 



l[\det{A, 



where Aj — {a'fj ) is the adjacency matrix of Dj each entry a'fj of which equals the 
sum of weights of all arcs from vertices vf^ to v^^"^ in Dj, and aj — ctj{G) satisfies 
(l). 



Corollary 3.2. Let G be a bulk plane bipartite lattice of order N with 2n-natation 
symmetry (where N ^ 00). If — is a constant, then the entropy of G 



lim — logAf(G) 



N't: 



log|detp)|d,, 



where D is the matrix of order N' obtained from Aj by replacing each aj in entries 
of Aj with cos x + i sin x, i"^ = —1. 

Proof. By Theorem 3.1, we have 

n-l 

M(G)=\{\dei{A.j)\, 



where Aj is the adjacency matrix of Z), .Hence 



lim — log M(G) = lim 



2nN' ^ 



^ n — 1 

log I det{A, ) I = hm — ^ log I det{A, ) | . 



n^oo nN' ' 
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By the definition of Aj's, it is not difficult to see that 



lim - V log I det(A,)| = - / log I det(D)|d, 
j=0 -^0 

which implies the corollary. □ 

Similarly, we have the following 
Corollary 3.3. Let G be a bulk plane bipartite lattice of order N with 2n-natation 

N_ 

2n 

\og\Aet{D)\ 



symmetry, where N' — oo. Then the entropy of G 



lim I- log M(G') = - 
where D is the same as in Corollary 3.5. 



lim 



N' 



dx, 



4. TWO TYPES OF TILINGS OF CYLINDERS 

Two bulk lattice graphs, denoted by Gl{m,2n) and G2(m, 2n), are illustrated 
in Figure 5(a) and Figure 5(b), respectively, where G*(m, 2n) is a finite subgraph 
of an edge-to-edge tilings of the plane with two types of vertices — 8.8.6 and 8.8.4 
vertices, and Gl{m, 2n) is a finite subgraphs of 8.8.4 tilings in the Euclidean plane 
which has been used to describe phase transitions in the layered hydrogen-bonded 
SnCf ■ 2H2O crystal [25] in physical systems [2 [21 US]. The bulk lattice graph 
G*(m, 2n) is composed of 2mn hexagons whose fundamental part is a hexagon. 
Similarly, The bulk lattice graph G2(to, 2n) is composed of 2mn quadrangles whose 
fundamental part is a quadrangle. Physicists call each of this kind of bulk graphs 
"an (m, 2n)-bipartite graphs with the free-boundary condition whose fundamental 
domain is G" (see [T5]V 



3j\ 

a,- 





Figure 5. (a) The lattice graph G*(m, 2r7,). (b) The lattice graph G2(to, 2n). 



If we add edges (a^, a*), (&i, 5*) for 1 < z < m and (ci,c*) for 1 < i < 2n in 
Gl{m, 2n), we obtain an (m, 2n)-bipartite graph with the doubly-periodic condition 
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on a torus (see the definition in [Mj), denoted by G\{m,2n). Similarly, if we 
add edges {ai,a*) for 1 < m and {hi,b*) for 1 < j < 2n in Gl{m,2n), then an 
(to, 2n)-bipartite graph with the doubly-periodic condition on a torus, denoted 
by G2(to, 2n), is obtained. For some related work on the plane bipartite graph 
with the doubly-periodic condition, see Kenyon, Okounkov, and Sheffield [M] and 
Cohn, Kenyon, and Propp p]. Salinas and Nagle .25j showed that the entropy of 
G|(to, 2n), denoted by lim log[M(G|(m, 2n))], equals 



1 



7r/2 



5 + V25- 16cos2 



de w 0.3770. (2) 



In this section, as applications of the product theorem in Section 3, we enumer- 
ate perfect matchings of two types of subgraphs of tilings of cylinders — Gi(m, 2n) 
and G2(to, 2ri), where Gi(m, 2n) (resp. G2(to, 2n)) is obtained from Gl{m,2n) 
(resp. G2(to, 2n)) by adding extra edges {ai, a*), {bi,b*) for 1 < i < to (resp. 
(fli, o*), {bj,b*) for 1 < i < TO, 1 < j < 2n) between each pair of opposite vertices 
of both sides of them. We call each of Gi(to, 2n) and G2{m,2n) " an (m, 2n)- 
bipartite graph with the cylinder-boundary condition" (simply cylinder). We also 
obtain the exact solutions for the entropies of two types of corresponding tilings 
of cylinders. We observe that both the Archimedean 8.8.4 tilings of the Euclidean 
plane and the corresponding tilings of the cylinder have the same entropy which 
is approximately 0.3770. Based on the result obtained in [14], we also show that 
both Gi(to, 2n) and G{{m, 2n) have the same entropy which approximately 0.3344. 
These are reasonable conclusions from the physical intuition. 

4.1. THE CYLINDER Gi(m,2n). Note that the cylinder Gi(m,2n) can be re- 
garded as a 2n-rotation symmetric plane bipartite graph. There exists a regular 
orientation Gi(to, 2n)'^ of Gi(m, 2n) as stated in Section 3, which is illustrated in 
Figure 6(a). For Gi(m, 2n)'^, all hexagons in the first column have the same ori- 
entation, all hexagons in the second column have the inverse of the orientation of 
hexagons in the first column, and so on. 

Theorem 4.1. For the cylinder Gi(to, 2n), the number of perfect matchings of 
Gi(to, 2n) can be expressed by 

1 1 \ , I ^ 2m+l , I ^ 2m+l' 

M(Gi(TO,2n)) = -n^— (V4^ + /5.) +(V4^-/5.) 
j=0 ^4 + L 

(3) 

and the entropy of Gi{m,2n), i.e., lim -p^ — logM(Gi(m, 2n)), equals 

m,n— *oo -'-^^^ 

2 

— / log(cosa; -I- V4 + cos^ a;)(ix « 0.3344, 
Stt Jo 

where j3j = cos ^ if n is odd and fij ~ cos ^^"'^^^^ otherwise. 

Proof. For the regular orientation Gi{m,2nY of the cylinder Gi(to, 2n) shown in 
Figure 6(a), by the definition of Dj for < j < n — 1 defined above, the correspond- 
ing digraphs -Dj's have the form illustrated in Figure 7(a), where a^'s satisfy (1), 
and each arc {vi,Vj) in Dj whose weight is neither aj nor must be regarded 
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(a) (b) 



Figure 6. (a) The regular orientation Gi{m,2nY of Gi(m, 2n). 
(b) The regular orientation G2{m,2ny of 6*2 (m, 2n). 




(a) (b) (c) (d) 

Figure 7. (a) The digraphs Dj's corresponding to the regular 
orientation Gi{m,2nY in Figure 6(a). (b) The digraph D'^ ob- 
tained from Dj by deleting vertex c{ shown in Figure 7(a). (c) The 
digraphs -Dj's corresponding to the regular orientation G2{m, 2nY 
in Figure 6(b). (d) The digraph D'^ obtained from Dj by deleting 
vertex h\ shown in Figure 7(c). 
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as two arcs {vi,Vj) and {vj,Vi) with weights 1 and —1, respectively. Let D'j be the 
digraph obtained from Dj by deleting vertex (see Figure 7(b)). 
For j = 0, 1, . . . , n — 1, set 

L„,(j) = det{A,), L'^U) = dct(4), 

where Aj (resp. A'j) is the adjacency matrix of the digraph Dj (resp. D'j). It 
is not difficult to prove that {Lm{j)}m>Q and {L'^{j)}m>o satisfy the following 
recurrences: 

LmU) = (4 + 4/32)L„_i(j) + 4(3,L'^_,ij) for m > 1, 



LoU) = 1, L'.ij) = 0. 
Hence we have the following: 



for m > 1, 



Lrnii) = (8 + 4/32)i:„,_i(j) - 16Lm-2{j) for m > 2, 



Lo(i) = 1, Li{j) 
which implies the following: 
1 



LmU) 



4 + /3?+A 



2m+l 



+ U/4 + /5 -A 



2m+l' 



(4) 



Hence the equality (3) follows from the product theorem (Theorem 3.1) and (4). 
So the entropy of Gi(m, 2n) 

2 1 
■logM(G'i(m, 2n)) = lim 



lim 

m,n — ^oo 



12m7i 



aog2-i^log(4 + /?|) + ^log (^4 



m.n^oo Qmn 
2m+l 



+ /3 +A 



lim tt'— log 



i=o 



m,n— J-CX5 



J-0 



(^IT^ + ft)"'"*' + (y^-,(J,) 

(v' 



lim y] log 



4 + /5f + A 



2m+l 



2m+l' 



2m+l 



lim 

m,n — 'oo 
2 f^/Z 



2m+l 
3mn 



J=0 



E log(J4 + /32+/3,) 



— ^ /J" log(cos a; + \/4 + cos^ x)dx 
and the theorem thus follows. 



0.3344 



□ 



In order to prove the following corollary, we need to introduce a formula of the 
entropy for an (n, n)-bipartite graphs with the doubly-period condition obtained by 
Kenyon, Okounkov, and Sheffield [14]. Let G be a Z'^-period bipartite graph which 
is embedded in the plane so that translations in the plane act by color-preserving 
isomorphisms of G — isomorphisms which map black vertices to black vertices and 
white to white. Let G„ be the quotient of G by the action of nZ^. Then G„ is a 
bipartite graph with the doubly-period condition. Let P{z, w) be the characteristic 
polynomial of G (see the definition in page 1029 in [11]). Authors in [14] showed 
that the entropy of G„ 



n-+oo n^lGil 



logAf(G„)= f f log|P(z,H|--, 



(5) 
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where D — {(z, w) £ : \z\ = \w\ = 1} and P = —1. 

Corollary 4.2. Both Gi(m, 2n) and G\(rn,2n) have the same entropy, that is, 
Urn — ^log(M(Gi(TO,2n))) Urn log(M(G*i (to, 2n))) w 0.3344. 

m.n^oo 12TOn m.n^oo 12mn 



Figure 8. The fundamental domain of G{{m,2n). 

Proof. Note that by the definition in [TJ] the fundamental domain of G\{m,2n) 
is composed of two hexagons (see Figure 8). Otherwise, if we use a hexagon as 
the fundamental domain, then it does not satisfy the condition "color-preserving 
isomorphisms" . It is not difficult to show that the characteristic polynomial of 
Gl{m,2n) 

P{z,w) = 10-4(z + z"i)- {w + w^^). 
Hence, by (5), we have 

hm ^_log(M(G*(™,2n))) 



_2 rSir r2Tr , 

(27r)2 Jo Jo 

24W lo'' l°g(10 - 8 COS x -2cosy)d^dy. 



12(1^)2 Jo'' Jo log(10 - 8 COS a; ~ 2 cos y)d^dy 

(•2-K [•2lT ■ 



Let F{y) = Jg'^log(10 — 8 cos a; — 2 cos y)dx — 2 JJ^log(10 — 8 cos a; — 2coay)dx 
Then 

dx 271 sin y 



F'{y) =Asiny [ 
Jo 



10 - 8cosx - 2cos?/ ^cos2 y - 10 cosy + 9 



Hence we have 



F{y) — 27r log(5 — cos y + \/ cos^ y — 10 cos ?/ + 9) 

implying 

-z--^ / log(10 - 8cosa; - 2cos?/)d2;dy 
247r^ 7o Jo 

= log(5 — cosx + V cos^ X — 10 cos x + 9)dx 

~ ^ lo ^°S(5 ~ ^ + V cos^ a:: — 10 cos x + 9)dx 
So, by Theorem 4.1, it suffices to prove the following: 

4 / log(cos X + \/4: + cos^ x)dx = / log(5 — cosx + \/ cos^ x — 10 cos a: + 9)dx. 
Jo Jo 

Note that 

4/ log(cos X + ^4 + cos^ x)dx = 2 / ln(2 cos^ x + 4+ \/l6 cos^ x + 4 cos'' x)dx 
Jo Jo 

= / log(2cos^ ^+41 /l6cos2 ^ + 4cos4 = / log(5+cosx+-\/cos2 x + 10 cos x + 9)dx 
Jo 2 V 2 2 Jq 
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= / log(5 — cos X + \/ cos^ X — 10 cos x + 'd)dx- 
Jo 

The corollary has thus been proved. □ 

4.2. THE CYLINDER G2(m, 2n). Note that the cylinder Gi{m, 2n) can be cm- 
bedded into the plane such that it is a 2n-rotation symmetric plane bipartite graph. 
There exists a regular orientation G2{m,2ny of G2{m,2n) as stated in Section 3, 
which is illustrated in Figure 6(b). 

Theorem 4.3. For the cylinder G2{m,2n), the number of perfect matchings of 
G2{m,2n) can be expressed by 

M(G2(m,2n)) = 



n 

j=0 



2 



9 + 16/3| V ^ / 2^9 + 16/3: 



(6) 



and the entropy, i.e., lim o:;^ logM(G2(m, 2n)), equals 

m.n— »(X) 



1 

2^ 



7r/2 

log 







5 + \/25- 16cos^6' 



« 0.3770, 



2 

w/iere /3j = cos ^ if n is odd and Pj = cos ^^^2n^^ otherwise. 

Proof. We can prove easily the statement in the theorem on the entropy from (5). 

Hence it suffices to prove that (5) holds. For the regular orientation G2{m, 2nY of 
the cylinder G2(m, 2n) shown in Figure 6(b), by the definition of Dj for < j < n—1 
defined above, the corresponding digraphs D^-'s have the form illustrated in Figure 
7(c), where a^'s satisfy (1), and each arc {vi.Vj) in Dj whose weight is neither 
Uj (or —aj) nor aj"^ (or — Q;~^) must be regarded as two arcs {vi,Vj) and {vj,Vi) 
with weights 1 and —1, respectively. Let D'j be the digraph obtained from Dj by 
deleting vertex bl (see Figure 7(d)). 
For j = 0, 1, . . . , n — 1, set 

PUl) = dct{A,), P;„(,?) = dct(4), 

where (resp. Aj) is the adjacency matrix of the digraph £)j (resp. Dj ) illustrated 
in Figure 6(c) (resp. Figure 6(d)). It is not difficult to prove that {Pm{j)}m>o and 
{Pm{j)}m>o satisfy the following recurrences: 

P2m+l{j) = 4P2™(j) + 2p,P/,Jj), P^„,+ ,{j) = -2P,P2m{j) " P^^ij) m > 0, 

P2m{3) = ^P2m-l{j) " 2PjP^^_,{j), P^^{j) = 2/3,P2„_i (j) - P^„_i(j) m > 1, 

Po(j) = 1, P^(i) = 0, Pi(j) = 4, Pi{j) = -20J. 
Hence from the recurrences above we have the following: 



Wi(i) \_( 16 + 4/3j -10f3j \ f P2m-i{3) 



P2m+l{j) J V -lOA- 4/32 + 1 ; V P^^_,(j) 

and 

P2m{j) \_f 16 + 4/32 10/3, \ / P2(„-l)(j) 

PLU) J ~ V m 4/3| + 1 J 1, p^(„_i)(j) 



(7) 
(8) 
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Let am — P2m+i{j) and bm ~ P2m(S) fo'' 3 ^ 0- Hence we have 

a™ = (I8/32 + 17)a„_i - 16(1 - /?|)2a„_2 for m > 2, 
6,„ = (I8/32 + 17)6„,_i - 16(1 - 0\fh,r^-2 for m > 2, (9) 



ao 4, ai = 64 + 36/3|, 5o 1, 61 = 16 + 4/3|. 
By solving the recurrences in (9), then we have 



and 



n-l 

n 

J=0 



ri-1 



y^9 + 16/3| + 3 /s + y^9 + 16/3|\ y^9 + 16/3^ - 3 /s - .^9 + 16/?: 



2^9 + 16/32 V ^ / 2^9 + 16/3^^ 



6m = n 



2m j . I . 2m' 



16/3^ + 3 ^ + Joi + I6/32 \ ^9 + I6/32 - 3 5-J9 + 16/3^ 



2^/9 + 16/32 V / 2^9 + 16/32 



Hence (6) has been proved and the theorem thus follows. □ 

From (2) and Theorem 4.2, both the 8.8.4 tilings of the Euclidean plane and the 
corresponding tilings of the cylinder have the same entropy which is approximately 
0.3770. 

5. CONCLUDING REMARKS 

In the classical work on the dimer problem for plane quadratic lattices, Kasteleyn 
[llj proved that the m x n quadratic lattice with the free-boundary condition, the 
m X n quadratic lattice with the doubly-period condition (i.e., the m x n torus), 
and the m x n quadratic lattice with the cylinder-boundary condition (i.e., the 
mxn cylinder) have the same entropy (= 2c/7r « 0.5831, where c is Catalan's con- 
stant). In this paper we have investigated the problem on enumeration of perfect 
matchings of plane bipartite graphs with 2n-rotation symmetry. We compute the 
entropy of a bulk plane bipartite lattice with 2n-notation symmetry. We obtain 
the explicit expressions for the numbers of perfect matchings and entropies for two 
types of tilings of (the surface of) cylinders and we showed that both the bipar- 
tite graph Gi{m,2n) with the cylinder-period condition and the bipartite graph 
G* (to, 2n) with the doubly-period condition have the same entropy and that both 
the bipartite graph G2(m, 2n) with the cylinder-period condition and the bipartite 
graph G\{m, 2n) with the doubly-period condition have the same entropy. A nat- 
ural problem is how to enumerate perfect matchings of the plane bipartite graph 
with (2n -t- l)-rotation symmetry. A more general problem is to enumerate per- 
fect matchings of the plane graph (not necessary to be bipartite) with n-rotation 
symmetry. 
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